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For the fermion field in the two-dimensional Gross-Neveu model, we introduce a flow 
equation that allows a simple 1/N expansion. By employing the 1/N expansion, we 
examine the validity of a universal formula for the energy-momentum tensor which 
is based on the small flow-time expansion. We confirm that the formula reproduces 
a correct normalization and the conservation law of the energy-momentum tensor by 
computing the translation Ward-Takahashi relation in the leading non-trivial order in 
the 1/N expansion. Also, we confirm that the expectation value at finite temperature 
correctly reproduces thermodynamic quantities. These observations support the validity 
of a similar construction of the energy-momentum tensor via the gradient/Wilson flow 
in lattice gauge theory. 
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1. Introduction 

It has been well recognized 0,8 that the construction of the energy—momentum tensor 
the Noether current associated with the translational invariance—is quite involved in lattice 
field theory. This is because lattice regularization explicitly breaks the translational invari¬ 
ance and the energy-momentum tensor is a composite operator containing local products 
of held variables. Because of radiative corrections, a naive discretization as it stands cannot 
reproduce a correct normalization and the conservation law of the energy-momentum tensor 
in the continuum limit. Recently, a completely new approach to this problem, on the basis of 


, has been proposed 0- 


the gradient/Wilson how [3 h 5| and the small how-time expansion 

8 in the context of lattice gauge theory. In that approach, especially that in Refs. |YJ, iaj, one 
constructs a “universal formula” for the energy-momentum tensor using a perturbative solu¬ 
tion of the gradient how. This construction relies on the UV hniteness of the gradient how 


in gauge theory [4|, |6|] such that renormalization of any composite operator of flowed helds 
is very simple. The universal formula is supposed to provide a regularization-independent 
expression for the energy-momentum tensor and thus is expected to be usable even with 
lattice regularization. 

The above approach is based on natural assumptions such as the existence of the energy- 
momentum tensor and the renormalizability of the gradient how in the non-perturbative 
level. Also, the formula in Ref. 0] has been numerically tested for quenched QCD at finite 
temperature [101, 111]. However, it still remains important to investigate the validity of the 
approach in various possible ways. In particular, it is of great interest whether and how the 
universal formulas in Refs. [7|, |9j, which are constructed by using perturbation theory, can 
capture non-perturbative low-energy physics or not. 

As shown in Ref. [12], the gradient how in the two-dimensional 0(N ) non-linear sigma 


model 13j] possesses a UV hniteness quite similar to that of four-dimensional gauge theory. 
By utili zing this UV hniteness, one can imitate the above construction of the universal 
formula 121 ]. For the two-dimensional 0(N ) non-linear sigma model, the 1/IV expansion is 
available and, to some extent, the gradient how can also be solved in the large N limit 14.1151]. 


In Ref. [l4], using this non-perturbative solution, the universal formula for the energy- 


momentum tensor has been analytically tested by computing the expectation value at finite 
temperature. The expectation value correctly reproduces thermodynamic quantities obtained 
by the conventional 1/IV expansion. This study demonstrates that the universal formula 
reproduces a correct normalization at least for those quantities. 

Another interesting issue is whether the conservation law (and more general Ward- 
Takahashi relations associated with the translational invariance) is correctly reproduced by 
the universal formula. This analysis for the two-dimensional 0(N ) non-linear sigma model 
has not been carried out, because in Ref. |lj] the gradient how was solved only in the leading 


order in the 1 /N expansion with which any correlation function is factorized into one-point 
functions @ 

In the present paper, with the above motivations, we consider a similar universal for¬ 
mula for the energy-momentum tensor in the two-dimensional Gross-Neveu model 16[. The 
point is that, in this non-gauge, unconstrained system, one can introduce a very simple how 


1 It might be possible to use the large N solution given in Ref. [l5| to investigate this issue. 
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equation that does not contain any interaction. Although this is not the gradient flow in 
the sense that the flow is defined with respect to the equation of motion of the original 
system, such a choice is perfectly legitimate from the perspective of the UV finiteness of 
the flow. Similar simplification has also been adopted for the flow of the fermion field in 
gauge theory [3], Because of this simplification in the flow equation, the conventional 1/IV 
expansion H y directly provides the solution of the flowed fields. We can then readily 
examine, in the leading non-trivial order in the 1/IV expansion, if the universal formula 
correctly reproduces the translation Ward-Takahashi relation 
This paper is organized as follows. In Sect. [2] we introduce a flow equation for the fermion 
field in the Gross-Neveu model. In Sect. [3] along the line of reasoning in Refs. 0, SQ, 
we construct a universal formula for the energy-momentum tensor in the present system. 
This construction itself is based on one-loop matching with the expression with dimensional 
regularization. In Sect. 01 we recapitulate the conventional 1/IV expansion of the present 
system. Section[5]is the main part of the paper and, in the leading non-trivial order in the 1/IV 
expansion, we examine if the universal formula correctly reproduces (some particular cases 
of) the translation Ward-Takahashi relation. Here, we observe that the universal formula 
precisely reproduces expected relations with the presence of the non-perturbative mass gap, 
although the construction of the universal formula itself uses one-loop perturbation theory. 
As another support for the universal formula, in Sect. 01 we compute the expectation value 
of the energy-momentum tensor defined by the universal formula at finite temperature 
as Ref. 


141 ]. It reproduces the correct results. The last section is devoted to conclusions. 


2. Flow equation in the Gross-Neveu model 


The Euclidean action of the Gross-Neveu model 


161 ] is given by 


S 


d D x 


(x)^i/ 1 (x) 



( 2 . 1 ) 


where D = 2 for our target theory and the fermion field has N components (* = 1,2,..., 
N) 0 In this system, we introduce a flow equation. That is, we introduce a fictitious time t 
and suppose that the fermion field evolves according to 


d t x(t, x ) = dpduxix), x(t = 0,x)= 'ip(x), (2.2) 

dtx(t, x) = d^xix), x(t = 0,x) = ip(x), (2.3) 

where the initial value for the evolution is given by the original fermion field which is the 
subject of the functional integral (with the distribution defined by Eq. (12. ID ). Note that 
Eqs. (12.21) and (12.31) are very simple; the right-hand sides are defined by the free Laplacian 
without any interaction. This is possible for the present non-gauge, unconstrained system. 
Although the above flow is not the gradient flow in the sense that the flow is defined by the 
equation of motion for the original action p.ll) . such a choice is completely legitimate as far 
as a UV finiteness of the flow is concerned—see the following discussions. 


2 The summation over repeated “flavor” indices i, j, ..., and Lorentz indices fi, v, ..., is always 
understood in this paper. 
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Since flow equations (12.2j) and f)2.3[i do not contain any interaction, the flowed fermion 
field becomes a simple linear functional of the fermion held at the zero flow time. That is, 


X(t, x) 


/ 


d D yK t (x -y)t/>(y), x(t,x) 


/ 


d D yK t (x-y)ijj(y), (2.4) 




is the heat kernel for the free Laplacian. Using Eq. (12.41) . correlation functions of the flowed 


fermion held can directly be obtained in terms of correlation functions of the original fermion 
held. For example, since the tree-level propagator of the original fermion held is given by 



(2.7) 


the tree-level propagator of the howed held is 



( 2 . 8 ) 


Also, the renormalization property of the unflowed fermion held is directly inherited by 
the howed fermion held. In particular, their wave function renormalization constants are 
identical. This is quite different from the howed fermion held in gauge theory B in which 
the wave function renormalization constant for the howed fermion held is independent of 
that of the original fermion held, due to interaction in the how equation. 

3. Universal formula for the energy-momentum tensor 

In this section, following the idea of Refs. B B0 , we construct a universal formula for 
the energy-momentum tensor in the Gross-Neveu model (12.11) using the small how-time 
expansion [b]. We hrst assume dimensional regularization with D = 2 — e and derive the 
explicit form of the energy-momentum tensor. Since dimensional regularization preserves the 
translational invariance, that energy-momentum tensor fulfills the Ward- Takahashi relation 
associated with the translational invariance; this implies that the energy-momentum tensor 
is correctly normalized and is conserved. However, the energy-momentum tensor with dimen¬ 
sional regularization is useful only in perturbation theory. Our universal formula below is 
intended to provide a regularization-independent expression for the energy-momentum ten¬ 
sor. This universal formula is thus also expected to be usable with lattice regularization for 
example, with which non-perturbative calculations are possible. 


3 Throughout this paper, we use the abbreviation 



(2.5) 
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Assuming dimensional regularization, the energy-momentum tensor can be obtained from 
the variation of the action 

6S = - J d D xUx)d»T^(x) (3.1) 

under the transformations 


H(x) = ^{x)8^{x), 5ip(x) = ^{x)d^{x). 
The explicit form is given 


(3.2) 


1 . 


T^(x) = ^ l {x) + ip l (x) - 5^ <j ^ l (x)^^ l (x) - ^- [i/j z (x)ip l (x)] 




1 ,,, 


Aq 


2 N 

(3.3) 

^_ y ^_ ' 

where 8^ = 8^ — 8 This operator does not receive the multiplicative renormalization, 
because of the translation Ward-Takahashi relation 


^C>ext jf d D X 8 ll T llu (x) O mt ^ = - (Oex 1 8 u O- mt ), 


(3.4) 


where V is a bounded integration region, C> ext is an operator outside the region and 0- mt 
is an operator inside the region. We define a renormalized energy-momentum tensor by 
subtracting the (potentially UV-divergent) vacuum expectation value as 


(x) — T /lu (x) (T^v^x)) . 


(3.5) 


Now, to derive the universal formula, we express the composite operator m in terms 
of the composite operator of the flowed fermion field. This can be archived by the so-called 


small flow-time expansion in Ref. [6| . By a one-loop perturbative calculation similar to that 
of Refs. 


, 12j], we find 


«->- 




X l (t,x) 8 v + Xv 8 Ax\t,x) - ( x\t,x ) 7 m d u + lu 8 Ax\t,x) 




4-7 


= i/>'(x) + xip l {x) - 


+ ln(87rt) + 1 


5^[^\x)il:\x)\ +0(f) (3.6) 


and 


[x\t,x)x'(t,x)] 2 - l[x'{t,x)x'(t,x)} 2 


7r 


—b ln(87rt) 


[^\x)^\x)] 2 + 0(t). 


(3.7) 


In this and following one-loop computations, we retain only terms leading in the large N 
limit, because only leading terms are relevant in the analyses in the following sections. 


4 Here, we have taken only the part of the expression appearing in Eq. being symmetric 

under /ifH'; the anti-symmetric part generates the Lorentz transformation and is not explicitly 
considered in what follows. 
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From Eq. (|3.6p . we also have 


X z (t,x)^x l {t,x) - (x\t,x)*$x\t,x) 


= ^*(x)VV(*) - 77 — 

iV 7T 


ln(87rt) 


[i/; l {x)ip l (x)Y + 0(t). 


(3.8) 


The relations (13.61) . (13.71) . and (13.81) may be inverted for composite operators of the unflowed 
fermion field. We then substitute those expressions in Eq. (|3.3I) to yield 


4"A 


4"A 


{ t iiAr ( x ) = 7 x\t, x) d v + Xv d A x*(t, x) - -S^x^t, x) ft x*(t, x) 


,4-A 


Ao 


2 N 


Ao 


+ 7777 1 + 7^ 


27r 


—b ln(87rt) + 1 


<V x )x i (ti x )] 2 — VEV 


+ o(t), 


(3.9) 


where VEV denotes the vacuum expectation value of the composite operator appearing in 
the right-hand side. In the one-loop order, the coupling constant is renormalized in the 
minimal subtraction (MS) scheme as 


Aq — //A ( 1- 


7T e 


(3.10) 


Then in terms of the renormalized coupling A, we have 


1 




<bA 


{T^}r( x ) = c iAiA~:X l (t, x ) (ill <9 ^ + lu d A x l (t,x ) - c 2 (A; /j)-(5 f x'(1, x) ft x*(t,x) 


. VA 


where 


+ c 3 (A; fj,)Sn V [x l (t, x)x l (t, x)] 2 - VEV 
+ 0(t), 

ci (A; n) = c 2 ( A;/x) = 1 + 0( A 2 ), 

cs(A; /i) = ^ 1 1 + ^ [ln(87 Tfi 2 t) + 1] | . 


(3.11) 

(3.12) 

(3.13) 


As we have noted, in the present system, the flowed fermion field receives the wave func¬ 
tion renormalization common to the unflowed fermion field. Since the fermion field does 
not receive the wave function renormalization to the one-loop order in the present system^ 
even composite operators of the bare flowed fermion field are UV finite without multiplica¬ 
tive renormalization; the flow ensures this UV finiteness. Then Eqs. (|3.12l) and (13.131) show 
that the right-hand side of Eq. (13.111) is UV finite. This should be so, because the energy- 
momentum tensor (after subtracting the vacuum expectation value) in the left-hand side 
must be UV finite. 


5 This persists also in the leading order 1/N expansion that is relevant to our analyses below. Thus, 
in this paper, we do not need to consider the wave function renormalization of the flowed fermion 
field. In gauge theory, on the other hand, renormalization of the fermion field has to be taken into 
account; see Ref. [9|. 
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Finally, we utilize a renormalization group argument. We apply the operation 

^)o (3 ' U) 

to both sides of Eq. (|3. 11 f) . where the subscript 0 implies that the bare quantities are 
kept fixed under the derivative. Since the energy-momentum tensor (13.311 is entirely given 
by bare quantities and the composite operators in the right-hand side of Eq. (13.111) are 
also bare, we infer that (/n9/<9/x)oCj(A; n) = 0 for i = 1, 2, and 3. These equations say that 
Cj(A;/r) = Ci(X(q)]q) for arbitrary q, where A (q) is the running coupling in the MS scheme 
with the renormalization scale q. Since the renormalization scale q in Ci(X(q); q) is arbitrary, 
we may take q = l/\/&i by using the flow time t. Then, since A(l/\/8t) —>• 0 for t —>• 0 by the 
asymptotic freedom, the above perturbative computation is justified for t 0. In this way, 
we arrive at 


{T tw \ 


R 


(x) = lim 

t-> 0 


T^{t,x) - 


(3.15) 


where 


1 


<-)• 




Tfiuit, x) = -x l {t,x) v+jvd M J x\t,x) - -5^x l {t,x) ft x\t,x) 




+ 


A(l/\/ 8 1 ) 


2N 


1 + 


A(l/V 8 1 ) 


2vr 


(In 7r + 1) 


<W [x\t,x)x l {t,x)\ 


(3.16) 


This is our universal formula for the energy-momentum tensor. This is universal in the sense 
that it does not refer to any specific regularization; the composite operator in the right-hand 
side is a renormalized quantity that must be independent of regularization as far as the 
parameters are properly renormalized. 

We stress that our computation which led to Eq. (13.161) is purely one-loop. Although 
we retained only large N leading terms in one-loop coefficients, no non-perturbative 1/N 
expansion is invoked at this stage. In particular, the fermion is treated as massless. We stress 
this point because the intention of the present paper is to see how the formula (13.161) that 
is obtained by one-loop perturbation theory can capture non-perturbative physics. More 
specifically, we want to see if the idea that coefficients in the universal formula can be 
determined by perturbation theory while low-energy non-perturbative physics is contained 
in matrix elements of composite operators works or not. This is the idea for the construction 
of the lattice energy—momentum tensor in Refs. [?], [9L, fl3 |. 

4. 1/N expansion in the Gross-Neveu model 

Now, for the analyses in subsequent sections, we briefly recapitulate the well-known non- 
perturbative solution in the present system (12.11) . an expansion in powers of 1/N [16, 17]. 

For a systematic 1/N expansion, it is convenient to introduce an auxiliary field cr(x) and 
rewrite the action m as 


S = d D x 


ip 1 {x)ftift' 1 (x) + a{x)ip l (x\ft l (x) + ~(j(x ) 2 

2Aq 


(4.1) 


Z = 


Ida(a 


exp 


d D x cr(x) z + N TV In [ft + a{x)\ 


(4.2) 


If we first integrate over the fermion field, the partition function becomes 

N 
2Xv 

Since the exponent is proportional to N in this expression, in the leading order of the 1/N 
expansion, the integral over the auxiliary held can be approximated by the value at the 

7 























saddle point. The saddle point is specified by the stationary condition for the exponent, i.e., 
by the gap equation, 

—cr = tr f — -. (4.3) 

Ao Jpip + <r 

The momentum integration in the right-hand side requires regularization. If we use 
dimensional regularization with D = 2 — e, we have 


1 1 

Ao * 7 vr 


- - - In 
e 2 


eV 2 M 


47T 




(4.4) 

(4.5) 


where 7 is the Euler constant. This tells us that, setting 

Ao = p e ^Z, 

the renormalization factor is given by 

Z~ 1 = l + -- (4.6) 

7 r e 

in the MS scheme. In terms of the renormalized coupling A, the saddle point is expressed as 

a 2 = 47 re _ 7 A 2 , A = /re”^. (4.7) 

As Eq. (14.11) shows, this saddle point provides a non-perturbative mass gap for the (originally 
massless) fermion. Corresponding to Eq. (|4.6p . the beta function is given by 


P = M 


d_ 

djl 


X=-eX- 


A 2 


7 r 


and thus the running coupling in the MS scheme is 

27T 27T 


A (g) = 


(4.8) 


(4.9) 


ln(A 2 /q 2 ) ln[e'i'cj 2 /( 47 rg 2 )] 

To obtain the next-to-leading order corrections in the 1/IV expansion, we have to consider 
the integration over the fluctuation around the saddle point in Eq. (14.21) . So we set 

cr(x) = a + 5cr(x). (4-10) 

The expansion of the exponent in Eq. m is then 
N 


2A 0 


J d D xcr(x) 2 + ./V Trln [(f) + a(x)] 
j d D xcr 2 + N J d D x tr J + a) 

/ dDxfe(l)2 

d°x / A°ySa(,)Sa(y) j tr / ~ f \ ~ ^ + ~ 


N 

2Aq 


N 

~2 


+ 0(6a 3 ). 

(4.11) 


There is no 0(5cr) term because a is the saddle point. After the momentum integration and 
the parameter renormalization (14.51) . we have 
N 


d D x cr{xY + N Tr In [p + a(x)] 


2A 0 

N 

47T 


d D xa z - 


N 

47T 


J d D x J d D y5a(x)5a(y) J e ^ x ~^B{p 2 , a 2 ) + 0(6a 3 ), (4.12) 
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where 


m 


B(p 2 ,a 2 ) 


jp 2 + 4 < 7 2 

p 2 


In 


/ Vp 2 + 4ct 2 + V/? \ 

\ VF 2 + 40-2 - V^V 


From this, the propagator of the fluctuating field 6cr(x) is given by 


(4.13) 


(Sa(x)Sa(y)) = '^ B(p 2 , a 2 )- 1 + 0(1/N 2 ). (4.14) 

In the above computation, we may use lattice regularization as well [3, Hi]- With the 
lattice spacing a, one may discretize the action (14.11) by replacing the Dirac operator fy by 
the Wilson Dirac operator for example, 


\ |>(^ + d l) - ad^d^] + m 0 , 


(4.15) 


where d^ and d* are forward and backward difference operators, respectively; mo is the bare 
mass parameter to be tuned to restore the chiral symmetry explicitly broken by the Wilson 
term. Then setting 

a = a + mo, (4-16) 


the gap equation with lattice regularization reads 


d 2 p 


a lat(* m » ) - tr 7 B (2»P^+Ia ! i2 + j 


= 0.7698-rin(a 2 <r 2 ) + 1.11861] 

a 2i 1 J 




(4.17) 


where Aq AT is the bare coupling with lattice regularization, B is the Brillouin zone B = { p^ 
—7T /a < p^ < ir/a] and 


_ 2 . (ap^\ 0 _ 1 . 

P,i = ~ sm J , Pf, = - sm (ap M 


(4.18) 


We choose the bare mass parameter mo so that the gap equation possesses a “symmetric 
solution” <j = 0; this corresponds to a massless fermion because a provides the fermion mass 
in the leading order of the 1/N expansion. This requirement leads to 


\ LAT 

m 0 = —0.7698 ——. 

a 

Under this choice, Eq. (|4.17|) says that 

cr 2 = e -l-1186 le -2VA& AT J_. 

a 1 

a has the same physical meaning as a in Eq. (14.71) . Thus, by choosing A^ at in Eq. (14.201) 
so that (7 = cr as a function of a, and rewriting everything in terms of this renormalized 
quantity, the dependence of physical quantities on adopted regularization disappears. In 
particular, it can be directly seen that the expression (14.121) also remains the same for 
lattice regularization (with a = cr). In what follows, we assume that this sort of parameter 
renormalization is made. 


(4.19) 

(4.20) 
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5. Restoration of the translation Ward Takahashi relation 

By using the large N solution in the previous section, we now consider correlation functions 
which contain the composite operator (13.161) . Then, by studying the small flow-time limit 
of the correlation functions, we examine if the energy-momentum tensor defined by our 
universal formula, Eq. (I3.15|) with Eq. (13.161) . fulfills (some particular cases of) the translation 
Ward-Takahashi relation, Eq. (|3.4I) . 

We first note that the fermion propagator in the 1/N expansion is, from Eqs. (14.11) 
and (14.101) . 

(^(x)^'(y)> = S ij [ j p( - x ~ y) + 0(1/N) (5.1) 

and thus the propagator of the flowed fermion field is given by 

( X i (t,x)x j (s,y)) = 8 i i [ e ip ^^^ + 0(l/N) (5.2) 

Jp Ip ~t (T 

by Eq. (12.41) . The propagator between the flowed and unflowed fermion fields, (%*(£, x)^ (y)) 
for example, is given by simply setting the corresponding flow time zero (s = 0 in this 
example) in Eq. (15.21) . 

The first correlation function we consider is 

(dfxT IJ , v {t,x)^' l (y)'4> l {z)^. (5.3) 

In the leading non-trivial order of the 1/N expansion, there are two types of connected 
diagrams which contribute to this correlation function; both are of 0(N). These two types 
of diagrams are depicted in Figs. |T] and [2], respectively. The contribution of type I diagrams 



Fig. 1 Type I diagrams which contribute to Eq. (15.31) . The solid line is the fermion 
propagator and the blob denotes the composite operator (|3.16l) . 


in Fig. [T] is, for small t, 
[dpf^{t,x)^ l {y)i)\z) 


,-tp 2 


e ip(y-x) e iq(x-z) N _ 

/ p J q ^ & 


x i(-p + <?) J ^ [7 pi(p + q)v + Xui{p + q)p\ - + 


+ ACVV50 

2i r 


1+ Ad/yM) + 

Z7T 




e ~tq 2 
ij + a' 


(5.4) 
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X 


X 


Fig. 2 Type II diagrams which contribute to Eq. (15.31) . The solid line is the fermion 
propagator and the blob denotes the composite operator ([3.161) . The broken line denotes 
the propagator of the auxiliary field, Eq. (14.141) . The interaction vertex between the fermion 
field and the auxiliary field (denoted by the small filled circle) can be read off from Eq. (|4.1I) 
with Eq. (14.101) . 


The second diagram in Fig. Q] contains a loop integral arising from the self-contraction in 
the last four-fermi term of Eq. (|3.16|) . The loop integral is finite, however, because of the 
Gaussian damping factor in the propagator ()5.2j) . We can rewrite the integrand in Eq. (|5.4I) 


as 


K~P + 7 [7 Ap + 9 ) 1 / + 7 'J(p + q)ii] - 77<W(^ + ! 


+ Mh§t ( 2eV t ) 

27T 


1 + S(i/yg) ■ 


2vr 




= Hi> + a) 


-iq u + + i),lv\ 


+ 


iPv ~ + i),lv\ 


Hi + a) 


+ ryi^) 

[ 27T 


1 + W^W+l) 


2vr 


+ 1 > i(-p + q) v a. 


On the other hand, from Eq. (|4.9j) . we have 

A(l/\/8 1) 


2vr 


and we find the following t —> 0 limits: 

t —^0 27 r 

Km hlZTM) 

t —>-0 27 T 

7 2 

p I rr 


ln(2e 7 cr 2 t/7r) 


ln(2e 7 <r 2 t) = — 1, 


lim ln(2e 7 cr 2 t) 
t —>0 


2e V t ) 


2vr 


= — In 7 r. 


lim 

t —>-0 


A(l/\/8 t) 


2i r 


ln(2e 7 cr 2 t) 


= 1. 


(5.5) 


(5.6) 


(5.7) 

(5.8) 

(5.9) 

(5.10) 


From these, we see that the last line of Eq. (|5.5j) vanishes for t —>• 0. Then when Eq. ([5.51) 
is substituted in Eq. (15.41) . the factor (i$ + a) in Eq. (15.51) cancels the external propagator 
1 /Hi) + (t) and then the integration over p produces the delta function S 2 (x — y) for t —> 0. 
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The situation is similar for the factor (ij + a) in (| 5.5 jl . In this way, we have 


Jim ^ d^Tuu, ( t , x)ip l (y)ip l {z)^ 

= -5 2 {x - y) (d u i/j l (y)ip l {z)) - 5 2 {x - z) {ip 1 {y)d v i> 1 {z)) 

+ O(N 0 ). 
(5.11) 

This is precisely the expected form of the Ward-Takahashi relation associated with the 
translational invariance. In fact, by considering the integration over the position x over the 
region that contains the points y and z, we observe that Eq. (13.41) with O m t = ip 1 (y)V’* (z) 
(and O ext = 1) holds @ 

Since the correct Ward-Takahashi relation is already saturated by type I diagrams, 
Eq. (j 5.11 j) . the type II diagrams in Fig. [2] should not contribute to the translation Ward- 
Takahashi identity. To see this, and for a later use, it is useful to compute first the left-hand 
side parts of the type II diagrams depicted in Fig. [3] An explicit computation of the diagrams 


+ d u 


& 2 {x - y)^hn,lv] {^ l {yW{z)) - 6 2 (x - z) {^\y)^\z)) ^[ 7,0 7*/] 



Fig. 3 The left-hand side parts of the type II diagrams in Fig. [2j 


in Fig. [3] yields 




- V + 2-£r 


-V(1 + ^M2 eVf) 


2vr 


1+ M±m } 

Z7T 


[ln(2e 7 a 2 i) + B(r 2 ,a 2 )] 


Using Eqs. (15.71) (15.101) . we then have 

lim Eq. flEED = ^ (<W ~ [B(r 2 ,a 2 ) - 2] . 

From this, for the type II diagrams in Fig. [2l 


(5.12) 

(5.13) 


t->o 


! in7 ( T iiv{ t i x W(y)' i l jl {z) 


ii 


P 


My-x) e iq(x-z) N(J ^ ^ _ 2B(r 2 ,a 2 )~ 1 ] - 


1 1 


if + g i<h + cr ’ 


(5.14) 


6 The energy-momentum tensor always possesses the ambiguity that results in the total divergence 

in the (unintegrat ed) W ard-Takahashi relation associated with the translational invariance. The 
second line of Eq. (15.111) , which corresponds to the Lorentz rotation generated by the anti-symmetric 
part of the canonical energy-momentum tensor, being the total dive rgen ce, does not contribute to 
the integrated form of the translation Ward-Takahashi relation, Eq. (1741) . 
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where 


r = — p + g, 


(5.15) 


and we have the desired result 

Jim =0. (5.16) 

Thus, in the leading non-trivial order of the 1/fV expansion, we have confirmed that the 
universal formula, Eq. (13.151) with Eq. (13.161) . reproduces the translation Ward Takahashi 
relation (13.41) for the product of elementary Helds, 0- mt = ifi 1 (y)'ip' 1 (z) (and O ex t = 1). This 
shows that the universal formula reproduces the correct normalization and the conservation 
law for the energy-momentum tensor, at least in the correlation function with elementary 
fields. 

The above computation in fact demonstrates that Eq. (13.41) is also reproduced for the 
scalar density operator, that is, 

Oint = ZsftivWiy), ©ext = 1, (5.17) 

where Z$ is an appropriate renormalization factor for the scalar density. In the leading 
non-trivial order of the 1/N expansion, there are two types of diagrams which contribute to 

, (5.18) 

as depicted in Figs. [I]and[5l 



Fig. 4 Type I diagrams which contribute to Eq. (|5.181) . In each diagram, the left blob 
denotes the composite operator (13.161) and the right blob denotes the scalar density operator 
in Eq. (15.171) . 



Fig. 5 Type II diagrams which contribute to Eq. (|5.18D . In each diagram, the left blob 
denotes the composite operator (|3.16D and the right blob denotes the scalar density operator 
in Eq. (15.171) . 


For type I diagrams, the computation is identical to that for Eqs. (15.121) and (| 5.13 [) . As 
is clear from Eq. (15.131) . we have lim^o^T^f, x)^ 1 {y)^ 1 {y))\ = 0. For type II diagrams 
also, we do not need a new calculation because Eqs. (15.121) and (|5. 131) (which correspond 
to Fig. m give the parts of the diagrams in Fig. 0 Thus, again from Eq. (15.131) . we 
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have lim.t^o{dfjTfj U (t, x)^ 1 (y))n = 0. These reproduce Eq. (13.41) with Eq. (15.171) because 

{d v \i(y)i/j l (y)]} = 0 by the translational invariance. 

We can further argue that Eq. (|3.4I) is reproduced when O ext is a collection of renormalized 
composite operators of the fermion field and 0 \ n t = 1. That is, for this situation, we can argue 
that 

Jim ^e>ext j d D x d/iTuvit, x)\ = 0. (5.19) 

This shows that the conservation law of the energy-momentum tensor is reproduced in 
the correlation function with generic composite operators. The argument is simple: There 
exist two types of diagrams which contribute to Eq. (|5.19|) . For type I diagrams in Fig. [6l 
we can use the identity (15.51) for fermion lines starting from the vertex of the composite 
operator (13.161) . Then, as Eq. (15.111) . we have 

Jim (y\yW{z) ■ ■ ■ dfxTfxuit, 

= -5 2 {x - y) (d^iyWiz) • • • ) - S 2 (x - z) (y)d v ft (z) ■■■) - 

+ [ 5 2 {x — y), S 2 (x — z), ..., inside the total divergence in x] . (5.20) 

Then for i / i / z, ..., the right-hand side vanishes. For the type II diagrams in Fig. [71 
from Eq. (j5.13|) . we simply have lim t _ > .o(C , ext^7) tl/ (t,x))n = 0. These imply Eq. (15.191) . 




Fig. 6 Type I diagrams which contribute to Eq. (15.191) . In each diagram, the leftmost 
blob denotes the composite operator (13.161) and other blobs denote the fermion composite 
operators contained in C> ext in Eq. (15.191) . 



Fig. 7 Type II diagrams which contribute to Eq. (15.191) . In each diagram, the leftmost 
blob denotes the composite operator (13.161) and other blobs denote the fermion composite 
operators contained in O ex t in Eq. (15.191) . 
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6. Expectation value at finite temperature 

The Ward Takahashi relation (15.111) shows that our universal formula for the energy- 
momentum tensor gives rise to the correct normalization at least within the correlation 
function with elementary fields. To give a further support on the correct normalization, in 
this section we compute the expectation value of the composite operator (13.161) at finite 
temperature and compare it with thermodynamic quantities directly obtained in the con¬ 
ventional 1/N expansion. A similar analysis for the two-dimensional 0(N) non-linear sigma 


model has been carried out in Ref. 


m- 


At finite temperature with inverse temperature f3, the propagator is given by 

e -(t+s)K+ P ?) 


^ n=—oo ** 

where u n is the Matsubara frequency 


lp(x-y) _ 


ijOUn + ijlPl + 07 3 


2irn 


bJn. — 


P 


+ 0(1/N), (6.1) 


( 6 . 2 ) 


and ap is the large N saddle point at finite temperature; ap is given by a finite-temperature 
counterpart of the gap equation (14.31) : 


J__ 1 f [dpi_ 

Ao J 2vr i 


1 


ilo^n + i'yiPi + crp 


(6.3) 


The required computation is almost the same as that in Ref. [141 ). because of the similarity 
of expressions. Using Eq. (|3.15l) with Eq. (13.161) . for the energy density e we have 


e = - ({T 00 } r (x))p 

N , n o. N 


47T 


(J1 2\ 
K (Tp - a ) — 


7 r 


a 2 pY2K 2 (/3apn) + O(N 0 


(6.4) 


n =1 


and, for the pressure P, 

P=({Tn} R (x))p 

n ( 2 2 , N 


^]T^ 2 (/%n) + O(A 0 


(6.5) 


71—1 


These are the results of the universal formula. 

On the other hand, the free energy density of the present system is given by 


N 


1 


f(P) = ivJ a p ~ ^ I ln K++ °f> ’ 


2A 0 


P 


2vr 


( 6 . 6 ) 


and the energy-density and the pressure are given by e = df(/3)/dB and P = —/(/?)//?, 
respectively. From comparison of Eq. (16.61) with Eq. (Al) of Ref. 1J], we see that these 


quantities can be obtained by making the substitutions /( 


-f(P), N ->• 2 N, A 0 ->• 2A 0 , 


crp —> ct|, and a —>• a 2 in Eqs. (All) and (A9) of Ref. [14(. We then observe a complete 
agreement with Eqs. (16.41) and (16.51) . This result again supports the validity of our universal 
formula for the energy-momentum tensor. 
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7. Conclusion 


The flow in quantum field theory and the small flow-time expansion can give rise to a 
regularization-independent expression for composite operators. In this paper, we examined 
the validity of a universal formula for the energy-momentum tensor by using the Gross- 
Neveu model and the non-perturbative 1/N expansion. In the leading non-trivial order in 
the 1/N expansion, we have observed that (some particular cases of) the Ward-Takahashi 
relation associated with the translational invariance is correctly reproduced by the universal 
formula even with the non-perturbative mass gap. This is interesting because the construc¬ 
tion of the universal formula itself requires only (one-loop) perturbation theory. We have 
also observed that the formula reproduces thermodynamic quantities correctly. These obser¬ 
vations support the validity of a similar construction of the energy-momentum tensor via 
the gradient/Wilson flow in lattice gauge theory. 
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Note added 


I learned that a strikingly analogous idea to define the energy-momentum tensor to ours 
has been considered on the basis of the operator product expansion (instead of the small 
flow-time expansion) in Ref. 201 ]. I would like to thank Jan Holland for the information. 
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